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ABSTRACT 

With observations of the EPCru system, we continue our series of measurements of spin-orbit angles in 
eclipsing binary star systems, the BANANA project (Binaries Are Not Always Neatly Aligned). We find a 
close alignment between the sky projections of the rotational and orbital angular momentum vectors for both 
stars (/3 p = -1.8 ± 1.6° and |/3 S | < 17°). We also derive precise absolute dimensions and stellar ages for this 
system. The EP Cm and DI Her systems provide an interesting comparison: they have similar stellar types 
and orbital properties, but DIHer is younger and has major spin-orbit misalignments, raising the question of 
whether EP Cm also had a large misalignment at an earlier phase of evolution. We show that tidal dissipation 
is an unlikely explanation for the good alignment observed today, because realignment happens on the same 
timescale as spin-orbit synchronization, and the stars in EP Cm are far from syncrhonization (they are spinning 
9 times too quickly). Therefore it seems that some binaries form with aligned axes, while other superficially 
similar binaries are formed with misaligned axes. 

Subject headings: stars: kinematics and dynamics -stars: early-type - stars: rotation - stars: formation - bina- 
ries: eclipsing - stars: individual (EPCrucis) - stars: individual (DIHerculis) - techniques: 
spectroscopic 



1. INTRODUCTION 

One might expect star-planet and close binary star systems 
to have well-aligned orbital and rotational angular momenta, 
since they originate from the same portion of a molecular 
cloud. However, there are also reasons to expect misaligned 
systems. Star formation is a chaotic process, with accretion 
from different directions at different times possibly leading to 
misalignment between the stellar a nd orbit rotation axes (e.g. 
Bate et al . |20 1 01 [Thies et al .1201 1 ) . There are also processes 
that could alter the stellar and orbital spin directions after their 
formation. For example a third body orbiting a close pair on 
a highly inclined orbit can introduce large oscillations i n the 
orbita l inclination and eccentricity of the close pair (Kozai 
|1962| ), thereby introducing large angles between the stellar 
spins and orbital angular momentum of the close pair. Close 
encounters and possible exchange of members in binary sys- 
tems (e.g. |Gualandris et al.||2004| ) would leave, among other 
clues, a fingerprint in the form of misalignment between the 
components. Tidal forces will over time erase these clues, be- 
cause dissipation will tend to bring the axes into alignment 
while also synchronizing the rotational and orbital periods 
(e.g. IZahn|| 19771 [Hut|p^T| [Eggleton & Kiseleva-Eggleton] 
|200l| ). Thus the degree of alignment between the stellar rota- 
tion axes and the orbital axis depends on its particular history 
of formation and evolution. Therefore measurements of stel- 
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lar obliquities - the angle between stellar equator and orbital 
plane - allow us to test theories of formation and evolution in 
close star-planet and star-star systems. 

For example the formation of star-star systems with orbital 
distances of only a few stellar radii is not completely under- 
stood. It seems unlikely that the stars formed at these or- 
bital distances because they would have overlapped during 
their pre-main sequence phase, when they had larger sizes. 
Therefore the orbital distance likely decreased after forma- 
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96% of binary stars with orbital periods less than 3 days have 
a third companion on a wide orbit while only 34% of bina- 
ries with orbital periods larger than 12 days have a third com- 
panion. Additional evidence for KCTF would be a misalign- 
ment between the stellar spin axes and the orbital spin, as- 
suming that close binaries have aligned axes at birth, and that 
tides have not had enough time to align the spin axes. Thus 
measurements of stellar obliquities in close binary systems to- 
gether with a good understanding of tidal dissipation in these 
systems might lead to a better understanding of binary forma- 
tion. 

For the case of close star-planet (hot- Jupiter) systems such 
an approach has been fruitful. Hot-Jupiters are thought to 
have formed much further from the star than their current or- 
bital distances, mainly because not enough material would 
have been available so close to the star. Different processes 
which could have transported the planet inward would lead 
to different spin-orbit angles, and indeed systems with both 
small and large spin-orbit angles have been found (see, e.g., 
Winn et al. |[2005l |Hebrard et al.||2008l |Johnson et al.| |2009 ; 
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TABLE l 
General data on EP Cru 



TABLE 2 
Observation Log for EP Cru 



HD 


109724 




NSV 


5783 




R-A.J2000 


12 h 37 m 17 s 


t 


Dec.nooo 


-56°47 / 17 // 


t 


Distance 


1.0(1) kpc 




Vmax 


8.7 mag 




Sp. Type 


B5V+ B5V 




Orbital period 


11?08 




Eccentricity 


0.19 




*P 


3.6(3)/?© 




Rs 


3.5(2)*© 




Mp 


5(1)M© 




M s 


4.7(1. 1)M 






15700(500)K 




T e ff s 


15400(500)K 





fData from ESA ( 1997} 

★Data from Clausen et al. 1J20071 

Note. — R p denotes the radius of the primary component and R s the 
radius of the secondary component. M p and M p denote the masses T e ff s and 
r e ff s denote the effective temperatures. 



ties. This suggests that the inward migration of hot Jupiters 
involves changes of the orbital planes of the planets. 

With the BANANA project (Binaries Are Not Always 
Neatly Aligned) we aim to get a better understanding of the 
formation of close binaries as well as their tidal spin evolu- 
tion. Here we study the EPCru binary system. This is the 
fourth system which we study as p art of the BANANA project 
dAlbrecht et al.||20Q7l|2Q09l|20TT| Pa pers I-III). We also re- 
fer the reader to Triauc TeT al.| (|2012 ), for a description of a 
similar project by other investigators. While most of the stars 
in our sample are of early spectral types, their EBLM project 
focuses on eclipsing systems harboring low-m ass stars. 

EPC ru was only recently characterized by Cla usen et al. 
(2007 ). Table [T] gi ves some syste m parameters. We selected 
this system because Clausen et al. ( 2007 ) found it to be similar 
to DI Her, for which we already foun d the spin-orbit angles to 
be very large ( Albrecht et al. 2009 ). In particular the orbital 
parameters, the stellar masses, and the projected stellar rota- 
tion speeds (vsin/*), are similar in these two systems. Here v 
indicates the equatorial rotation speed and i* the inclination of 
the stellar rotation axis along the line of sight (LOS). There is 
one important difference bet ween the two systems: the age of 
the stars. |Clausen et al.| ( |2007] ) estimated an age of « 50 Myr 
for the two stars in the EP Cru system while DI Her is essen- 
tially a Zero Age Main Sequen ce (ZAMS) syste m with an 
estimated age of 4.5 ± 2.5 Myr ( [Claret et al.|2010| ). Therefore 
by studying EPCru we have the opportunity to learn if ac- 
cording to our current understanding of binary evolution one 
system is simply an older version of the other, or if they had 
different childhoods altogether. 

The plan of this paper is as follows. In the following sec- 
tion we describe our observations. The analysis of the spec- 
troscopic observations during eclipses and out of eclipses is 
presented in Section [3] We present the results on the absolute 
dimensions, the orientations of the stellar rotation axes, and 
the derived age in section [4] For the remainder of the paper 
we focus on the interpretation of the obliquity measurements 
before we summarize our findings in the conclusions. 



Obs. mid time 
(BJDtdb) 


Phase 


Eclipse 


SNR 


2455192.85669 


0.70 




142 


2455193.83164 


0.78 




97 


A C C \ C\ A /I C\C f\C\ 

2455194.69509 


0.86 


sec 


77 


2455194.83360 


0.87 




107 


2455196.85095 


0.06 




107 


2455197.86201 


0.15 




115 


1HCC1AA crvT c c 

2455309.59755 


0.24 


pri 


84 


2455309.63874 


0.24 


pri 


86 


2455309.68141 


0.24 


pri 


91 


1 JCC1AA al OA/1 

2455309.71894 


0.25 


pri 


85 


2455309.75526 


0.25 


pri 


89 


1HCC1AA 

2455309.80483 


0.25 


pri 


81 


1HCC1AA O/IIOI 

2455309.84282 


0.26 




71 


2455309.86946 


0.26 




83 


2455310.46906 


0.31 




73 


A C C n 1 1 /I /'OTA 

2455311.46870 


0.40 




96 


2455313.48255 


0.59 




94 


2455314.57372 


0.68 




67 


2455338.45972 


0.84 


sec 


60 


i^ccno ,10011 

2455338.48823 


0.84 


sec 


82 


2455338.52361 


0.85 


sec 


94 


2455338.57536 


0.85 


sec 


97 


2455338.61565 


0.85 


sec 


62 


2455338.65529 


0.86 


sec 


64 


2455338.69515 


0.86 


sec 


103 


2455338.72942 


0.86 


sec 


54 


2455340.44257 


0.02 




78 


2455341.55552 


0.12 




140 


2455342.55950 


0.21 


pri 


120 


2455353.46096 


0.19 




130 


lyiCCO CO rAAHA 

2455353.50079 


0.20 




114 


1/ICCIf 1 

2455353.52531 


0.20 




144 


2455353.57109 


0.20 


pri 


123 


2455353.61561 


0.21 


pri 


137 


1/ICCIf 1 f C f\C f\ 

2455353.65050 


0.21 


pri 


128 


1/ICC1/A A n A c 

2455360.46345 


0.83 


sec 


83 


i/icc^tn coin/; 
2455360.53196 


0.83 


sec 


89 


2455360.60002 


0.84 


sec 


53 


2455360.64101 


0.84 


sec 


44 


2455361.49041 


0.92 




103 


2455362.49566 


0.01 




88 


2455628.85227 


0.06 




85 


2455630.66346 


0.22 


pri 


82 


2455630.72346 


0.22 


pri 


92 


2455630.78472 


0.23 


pri 


80 


2455630.85156 


0.24 


pri 


91 


2455632.71659 


0.40 




83 


2455635.68609 


0.67 




81 



Note. — The phase is defined such that phase corresponds to peri- 
astron. In column 3 'pri' indicates that the observation was taken during a 
primary eclipse and 'sec' indicates that the observation was obtained during 
an ongoing secondary eclipse. 



2. SPECTROSCOPIC OBSERVATIONS 

We observed the EP Cru system with the FEROS spectro- 
graph ( [Kaufer et aL]|1999| ) on the 2.2 m telescope at ESO's 
LaSilla observatory. We obtained 48 observations on mul- 
tiple nights between December 2009 and March 2011. Ta- 
ble[2]gives an observation log. The observations had a typical 
integration time of lOmin. On each night, ThAr exposures 
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Fig. 1 . — Spectra of EP Cru obtained at different orbital phases. Each panel shows spectra of both stars in the spectral region around the Mg n line. The 
orbital phase of the observation is indicated by the number in each panel, with phase occurring at periastron passage. The gray solid lines represent the data 
and the (red) dashed and (blue) dash-dotted lines are the simulated absorption lines of the primary and secondary, respectively. The black line is the best fitting 
model. The dots around the line at a flux level of 1 .05 represent the differences between the data and the model. 



were taken to calculate a wavelength solution and monitor any 
changes in the spectrograph. For all observations we used the 
MIDAS FEROS package installed on the observatory com- 
puters to reduce the raw 2D CCD images and to obtain stel- 
lar flux as a function of wavelength. The uncertainty in the 
wavelength solution, expressed in velocity, is few ms" 1 , and 
is negligible for our purposes. The resulting spectra have a 
resolution of w 50000 around 4481 A (the wavelength area of 
the spectra we analyze). We corrected for the radial- velocity 
(RV) of the observatory, performed initial flat fielding with the 
nightly blaze function, and flagged and omitted bad pixels. 

3. ANALYSIS 

In this section we outline the analysis of the spectra with the 
aim of deriving absolute dimensions of the system and learn- 
ing about the projected obliquities of both stars via measure- 
ments of the Rossiter-McLaughlin (RM) effect, which occurs 
during eclipses. We describe which part of spectrum we an- 
alyze and briefly introduce the model to which we compare 
the data and the algorithm used to extract system parameters. 
Our approach for EP Cru is similar to the approach employed 
in Papers I— III. 

Spectral region — We focus on the Mgll line at 4481 A, as 
this line is relatively deep and chiefly broadened by stellar ro- 
tation. It is located in the red wing of the pressure-broadened 
He I line at 4471 A. While this li ne might also be included in 
the analysis (Albrecht e t al.|20li) , we decided here to exclude 
it as there is enough signal in the Mg II line and modeling the 
pressure broadening in the He I line represents an additional 
complication. Thus we fitted a Lorentzian model to the en- 
croaching wing of the He I line and subtracted it before mod- 
eling the Mg II line. For this fit we used the spectral regions 
4472 A - 4476 A and 4486 A - 4498 A, thereby avoiding the 
influence of the Mg II line. Each spectrum was binned to a 



resolution of about 12 kins" 1 , to speed up subsequent com- 
putations. Because the stellar rotation speeds are an order of 
magnitude larger, there is no significant loss of information 
due to the binning; we verified this by experimenting with 
higher resolutions. 

Model — The measured spectra show absorption lines of both 
stars in the system. Before, during, and after eclipses the RVs 
of both stars are similar, leading to a substantial overlap of 
the two absorption lines. Hence, light emitted from both stars 
has to be accounted for when analyzing the RM effe ct. We 
used the numerical code from Albr echt et al.| ( |2007| ) which 
simulates the spectra of both stars in a system. 

The stellar disks are discretized with ~ 30000 pixels in 
a Cartesian coordinate system. We assume the stars to be 
spherical because they are well separated with rotation speeds 



much slower than the breakup velocity; Claus enet aL] ( [2007 ) 
estimates an oblateness of about 0.0008" We further as- 
sume uniform rotation and quadratic limb darkening. Stellar 
surface velocity field s are param etrized adopting the macro- 
turbulence model by |Gray| ( |2005| )p1 The coordinates of both 
stars projected on the sky are calculated and light from visible 
parts of the stellar hemispheres is integrated. The resulting 
absorption line kernels are shifted in wavelength correspond- 
ing to the line they represent and weighted according to the 
light contribution of the respective star. 

Parameter choices — Having the model in place we can now 
learn about the EP Cru system by specifying a number of pa- 
rameters. The Keplerian orbit of the two stars can be de- 
scribed with the following 6 parameters: The orbital period 
CP), a specific time of minimum light during primary eclipse 

7 Here we do not need to take the point spread function (PSF) of the spec- 
trograph into account as our binning (12 kms -1 ) is larger than the width of 
the PSF (6 km s- 1 ). 



4 



Albrecht et al. 




Fig. 2. — Spectra of EPCru obtained during primary eclipse. Similar to Figure[T] but this time for observations obtained during primary eclipses. The 
numbers on the left side of each panel indicate the observation mid-exposure times relative to the time of minimum light, in hours. Each inset shows an illustration 
of the eclipse phase of the background star. 




Fig. 3. — Spectra of EP Cru obtained during secondary eclipse. Similar to Figure[2] but for spectra obtained during secondary eclipses. 
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TABLE 3 
Results for the EP Cru system. 



Parameter 


This work 


Literature values 


Orbital parameters 


Time of primary minimum, T^ j (B JD-2 400 000) 


46181.6928zb0.0038 


46181.7068zb0.0003 1 


Period, P (days) 


11.0774707zb0.0000043 


11.0774701zb0.0000042 lc 


V^cosa; 


-0.104=b0.006 






0.4202zb0.0011 




ecosuj 


0.18187zb0.00034 


0.18162zb0.00007 1 


esinu 


-0.0450zb0.0026 


-0.0475zb0.0020 1 


Eccentricity, e 


0.1874=b0.0005 


0.1877zb0.0005 lc 


Argument of periastron, uo (deg) 


346.1zb0.7 


345.4=b0.6 1 


Cosine orbital inclination, cos i 


0.0052=b0.0014 




Orbital inclination, i Q (deg) 


89.70±0.08 


86.97zb0.09 lc 


Velocity semi-amplitude primary, K v (kms -1 ) 


102.2±1.5 


100.9il.3il0 1 


Velocity semi-amplitude secondary, K s (km s -1 ) 


106.2±1.4 


io5.9i3.5iio 1 


Velocity offset, j p (kms -1 ) 


-26.3±0.6 


-33 1 


Velocity offset, j p (kms -1 ) 


-27.2±0.5 


-33 1 


Orbital semi-major axis, a (Rq) 


44.83±0.37 


44.6zb4.2 1 


Stellar parameters 


Light ratio, L s /L p @ 4480 A 


0.8972±0.0020 


0.8972=b0.0020 lc 


Fractional radius primary, r p 


0.0801±0.0005 


0.0810zb0.0006 lc 


Fractional radius secondary, r s 


0.0779±0.0004 


0.0782zb0.0006 lc 


ul{+u2[ 


0.50±0.05 


(0.4+0.1)zb0.1 2c 


Macro -turbulence parameter, £ p (km s -1 ) 


22.3±1.7 




Projected rotation speed primary, vsim'p (kms -1 ) 


141.4±1.2±5 


150 1 


Projected rotation speed secondary, vsin/ s (kms -1 ) 


137.8zbl.lzb5 


150 1 


Projected spin-orbit angle primary, (3 P (°) 


-1.8±1.6 




Projected spin-orbit angle secondary, /3 S (°) 


< 17 a 




Primary mass, M p (Mq) 


5.02±0.13 b 


4.95zbl.06 1 


Secondary mass, M s (Mq) 


4.83±0.13 b 


4.72=bl.03 1 


Primary radius, R p (R@) 


3.590±0.035 b 


3.61ZC0.25 1 


Secondary radius, R$ (Rq) 


3.495±0.034 b 


3.48zb0.24 ! 


Primary log g p (cgs) 


4.028±0.008 


4.02i0.11 1 


Secondary logg s (cgs) 


4.035±0.008 


4U3iO.ll 1 



4.1 



Notes — 

a We prefer this larger uncertainty interval over the statistical uncertainty of -13 i 4° (Section 
b A solar radius of 6.9566 • 10 8 m was used 
c Value was used as prior. 

References — 

(l) |Clausen et"al]|2007) (2) |Claret|{2000 
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Fig. 4. — Orbit of EP Cru. (a) The apparent radial velocities (RV) in the EP Cm system for both stars as a function of orbital phase, defined such that phase 
zero occurs at periastron passage. The (red) filled circles indicate measured RVs of the primary and the open (blue) circles RVs of the secondary. Due to the small 
amount of light received from the eclipsed star during mid-eclipses it was not possible to assign RVs to the background stars for observations obtained within 
1 hr of mid-eclipse. The RVs shown here were not used in the analysis; they are displayed only for illustrative purposes, (b) Zoom-in on orbital phases close to 
primary eclipse. Only RVs of the primary are shown, (c) The same as b, but this time for the secondary eclipse, (d) A pole on view of the orbit of both stars. The 
line of sight (LOS) towards earth is indicated. 



(Tmin,i), the orbital eccentricity (e), the argument of perias- 
tron (oj), and the velocity semi-amplitudes of the primary and 
secondary stars (K{). Here the subscript 'i' stands for ei- 
ther 'p' indicating the primary star or Y indicating the sec- 
ondary (slightly less massive) star. In addition, velocity off- 
sets (7i) are needed. For e and oj we use the stepping param- 
eters y^coscj and y^sincj, as they are less correlated than e 
and oj themselves. Th e results from the photometric study by 
|Clausen et al.| ( |2007] ) can be used to constrain some of these 
orbital parameters. However the photometry used by [Clausen] 
et al.| ( |2007| ) was gathered about 20 years ago and the system 



is expected to have an apsidal motion period of a few thousand 
years. The change in oj could be of order 1° over the last two 
decades. In addition the apsidal motion is not measured yet 
and we cannot calculate it from the known system parameters 
as it depen ds on the true stellar obliquity and n ot only the sky 
projection ( |Shakura| 19851 |Albrecht et al.|2009| ). Thus we do 
not use the photometric values on T m \ n i : a nd oj as priors and we 
only use the measurements of P and e by |Clausen et al.||20 07 ) 
as prior constraints. We revisit this subject in section [43] 

Additional parameters are needed to describe the projected 
equatorial rotation speeds (vsin/i), the Gaussian width of the 
macro-turbulence (£), and the parameters of greatest interest 
for this study, the sky -projected spin-orbit angles (ffi). The 
angle is defined according to the convention of |Hosokawa| 
(1953}. 

The photometric character of the eclipses are specified by 
another set of parameters: the light ratio between the two stars 
at the wavelength of interest (L s /L p at 4480 A), the quadratic 
limb darkening parameters (ul{ and u2{), the fractional radii 
of the stars (ri), and the orbital inclination (i ), for which we 
step in cos/ . For the fractional radi i and the orbital incli - 
nation we use prior information from Clause net al.| ( |2007| ). 
For L s /L p we use their results in the b band. To constrain the 
limb darkening parameters we used the 'jktldfjtool to query 
the ATLAS atmospheres ( |Claret|2000| ) and placed a Gaussian 



http : //www. astro. keele.ac. uk/ jkt/ codes / jkt Id. html 



prior on ul{ + u2i with a width of 0.1 and held the difference 
ul{ - u2{ fixed at the tabulated value. 

An additional parameter is needed for each star to describe 
the relative depth of the Mg II lines. The Mg II line consists 
of a doublet, given the close spacing (0.2 A) we model it here 
as single line. 

The two components in the EP Cru system are very similar 
to each other (see Table [TJ. We therefore decided to use the 
same limb darkening parameters and macro-turbulence veloc- 
ities, for both stars, thereby reducing the number of free pa- 
rameters to 20. Of these, 7 are further constrained by Gaus- 
sian priors as explained above. Table [3] summarizes all of the 
prior constraints. 

There is always a small residual uncertainty in the initial 
normalization of the spectra. To propagate this into the un- 
certainty intervals of the final parameters we added for each 
of the 48 observation 3 free parameters which describe a 
quadratic function used to normalize the continuum level. The 
values of the normalization parameters were optimized us- 
ing a separate 3 -parameter minimization for each observation, 
each time a set of system parameters is evaluated. This pro- 
cess is similar to the "Hyper plane Least Squares " method that 
was described and tested by Bak os et aL] ( |2010| ). 

Parameter estimation — A MCMC code was used to obtain un- 
certainty intervals. The chains consisted of 0.5 million calcu- 
lations of x 2 . The results reported below are the median val- 
ues of the posterior distribution and the uncertainties intervals 
are the values which exclude 15.85 % of the values at each 
side of the posterior and encompassing 68.3 % of all values. 

4. RESULTS 

The results for the model parameters are given in Table [3] 
Figure [T] shows the spectra in the vicinity of the Mg II line and 
the corresponding model for the out-of-eclipse observations. 
Figures [2] and[3] show the same for the spectra obtained during 
primary and secondary eclipses. The apparent radial veloci- 
ties in the EP Cru system are shown in Figure [4] as well as a 
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Fig. 5. — Projected obliquity of the secondary star. Upper panel: A 
histogram of the posterior of /3 S with logarithmic y-axis. In addition to the 
main peak around an angle of -13° there is a smaller peak at similar but 
positive values. Lower panel: A random sup-sample of the posterior in the 
ft - r p plane. The horizontal lines indicate the photometric prior and uncer- 
tainty interval for r p . For larger values of r p smaller absolute values of /3 S are 
found. 

pole-on view of the orbit. 

4.1. Stellar Rotation and Projected Obliquities 

The main result of our analysis is that the sky projections of 
the two stellar rotation axes /3 p = - 1 . 8 ± 1 . 6° and /3 S = - 1 3 ± 4° 
indicate close alignment between the stellar rotation axes and 
the orbital angular momentum. However while the value for 
/3 P is consistent with prefect alignment /3 S seems to indicate a 
small but significant misalignment. How robust is this find- 
ing of a small misalignment? We note that we have some- 
what lower SNR observations during the secondary eclipse 
than during the primary eclipse and fewer observations di- 
rectly before, during and after the eclipse (See Table [2| and 
Figures [2j [3] andH]). To test the robustness of the result we 
reran the MCMC chain with different model assumptions. For 
example we constrained the model more, by leaving j{ and 
the line depths of both stars tied to each other, or we left £i 
and the limb darkening parameters completely free. We also 
excluded some observations to test if a small number of obser- 
vations are having a disproportionate influence on the result. 
For all these runs we found a negative /3 S . The result closest 
to alignment was /3 S = -9 ± 5°. At the same time the result 
for /? p did not change by more than 0.4° during these tests. 
There are, however, two peculiarities about our result for /3 S . 
The posteriors for all the other parameters have only a sin- 
gle peak, while the posterior of /3 S has a small (two orders of 
magnitude lower) secondary peak at positive angles around 
/? s = 13° (Figure [5] upper panel). In addition there is a corre- 
lation between /3 S and r p (Figure [5] lower panel). r p is the only 
parameter for which we do find a more than 1-a displacement 
between the prior constraints and results (Table [3]). Taken the 
above mentioned points into consideration we are confident 
that |/3 S | < 17° but we cannot exclude a small projected obliq- 
uity for the secondary star, with the data at hand. 
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Fig. 6.— Stellar evol ution tracks for EP Cru from the Yonsei- Yale series 
( Yi et al. 2001 Demarque et al. 2004) compared with the measurements. The 
tracks are interpolated to the measured masses and a metallicity of [Fe/H] = 
+0.03 that best fits the temperatures. The uncertainty in the location of the 
tracks that comes from the mass error is indicated with an error bar for the 
primary, and is similar for the secondary. The age according to these models 
is57±5Myr. 



For the projected rotation speeds we find v sin i v = 141.4 ± 
1.2 kins" 1 and v sin 4 = 137.8 =b 1.1 krns" 1 . We consider the 
formal uncertainties in the vsin/i to be too low for the fol- 
lowing reasons: We tested different limb darkening laws and 
found for example that for a linear limb darkening law the best 
fitting vsin/ values are lower by about 4 km s -1 . Also we sus- 
pect that our particular choice of parametrization of the stellar 
surface velocity fields will influence the values we find for the 
projected rotation. We therefore estimate that a uncertainty of 
5 km s -1 is more realistic and also indicated that uncertainty 
in Table [3] As mentioned above, normalization for each ob- 
served spectrum is included in our routine, hence any uncer- 
tainty in normalization is already incorporated in the formal 
uncertainty. 

We note that the projected rotation speeds are similar to 
the average rot ation speed for B stars (vsin / =130 krns -1 ), 
as analyzed by |Abt et aL| ( |2002| ). However the stars might 
have undergone a change of v due to tidal interactions (Sec- 
tion [5]). Therefore we can not conclude from the similarity 
of the measured vsin/ to the expected v that sin/ is close to 
unity. Nevertheless it seems unlikely that the stars have large 
inclinations relative to the line of sight and at the same time 
their projected axes on the plane of the sky are both aligned. 
In what follows we assume that not only the projections of 
the rotation axes are small, but that the axes themselves are 
aligned too (sin / ^ 1). 

Concerning the values for the macro-turbulence, we expect 
that the value we find does not have a simple physical inter- 
pretation. This is because we assume equal brightness of rais- 
ing and falling material as well as equal surface coverage of 
movement tangential and radial to the stellar surface, both as- 
sumptions do not need to be fulfilled in reality. We did test if 
there is a strong dependence of the measured values for pro- 
jected obliquities on our adopted model for macro- turbulence, 
and found none. 
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TABLE 4 

Expected apsidal motion in EP Cru 



TABLE 5 
EP Cru and DI Her 



Expected apsidal motion 


(arcsec per cycle) 


^GR 


1.87±0.023 


^Tides 


0.76±0.11 


<^Rot 


6.17±1.00 


^Total 


8.8±1.1 



Note. — o>gr denotes apsidal motion due to General Relativity, abides 
due to tidal distortions, and o>R 0t due to rotational distortion. 



4.2. Absolute Dimensions and Age 
From the posterior of our MCMC chain we find K v = 
102.2 ±1.5 km s -1 and ^ s = 106.2 ±1.4 km s -1 in agreement 
with values from the literature (Table [3]). We also calculated 
the K{ values only using out of eclipse data, making them less 
dependent on any assumption included in our eclipse model. 
With £p = 101.9 ±1.5 km s -1 and K s = 106.2 ± 1.6 km s~ l we 
obtain consistent results 

With the new spectroscopic data we not only obtain precise 
mass estimates for both stars but also improve on the absolute 
radii. This is because the accurate scaled radii r,- obtained via 
photometry need to be multiplied by the absolute scale of the 
system, the semi-major axis (a). 

With the new values for the stellar masses, surface gravi- 
ties (logg/), and the e ffective temperatures T e m measured by 
(Clau sen et al.||2Q07] see also Table [1} we can estimate the 
stellar ages using stellar evolution m odels. Here we employ 



the Yonsei -Yale evolutionary tracks ( |Yi et al.|2001| Demarque 

etallicity and an age 



|et al.|2004| ). We find a good fit for solar metallicity 
of 57 ±5 Myr (Figure [6]). Another good check is the tempera- 
ture difference between the stars, since the difference is prob- 
ably better determined than the absolute temperatures. Indeed 
the temperature difference predicted by the models (i.e., the 
separation between the evolutionary tracks) is in go od agree 



ment with t he temperature difference measured by |Clausen 
|etaTl ( |2007) . 

4.3. Apsidal Motion 

Now that the stellar rotation is known we can calculate the 
expected apsidal motion in the EPCru sys tem. We use the 
apsidal motion constant log(^) = -2.3 from |Claret| ( |2004| ) for 
both stars in the system. We assign an uncertainty of 0.1 in 
log space to this constant. From the results in Table [4] we 
can see that we expect a shift of ^ 1.6±0.2° over the last 20 
years (which approximately have elapsed since the photomet- 
ric measurements). Most of this shift is expected because of 
deformation of the stars by their rotation. That we find a small 
increase in the argument of the periastron and a (« 22 min) 
earlier primary eclipse than expected from linear ephemeris 
seems to indicate apsidal motion of the order of magnitude as 
expected. However spectroscopic data is not very good at de- 
termining uj and we find it difficult to estimate the significance 
of the measured value for uj. Therefore to make a meaningful 
comparison between the measured and expected apsidal mo- 
tion, new photometric eclipse timings should be undertaken. 

9 John Sou th worth provided us with the 5 out of eclipse spectra used in the 
Clausen et al. ( 2007) study and we found that these are consistent with our 
data set. Because of the potential small change in the argument of periastron 
over the last 20 yr they have not been included in this study. 



Parameter 


EPCru 


DIHer 


Sp. Type 


r$5 V+ dj V 


T5 C"\ J i T) C\ T-k 

dj V+ 135 V 


r (days ) 


1 1 no 
1 1 .Uo 


1U.JJ 




U.IO/ H-ZCL/.L/L/L/J 


480-1-0 nfn* 




5.02±0.13 


5.17±0.11* 


M s (Mo) 


4.83±0.13 


4.52±0.07* 


R P (Rq) 


3.590±0.035 


2.681 ±0.046* 


Rs (Rq) 


3.495±0.034 


2.478±0.046* 


A>(°) 


-1.8±1.5 


72=b4t 


A (°) 


< 17 


-84±8+ 


vsinz'p (kms -1 ) 


141.4±5 


108±4t 


vsinz' s (kms -1 ) 


137.8±5 


116±4t 


vsyn p (kms -1 ) 


16.40±0.16 


12.85±0.24 


vsyn s (kms -1 ) 


16.04±0.15 


11.89±0.24 


vps p (kms -1 ) 


19.51±0.19 


34.3±0.7 


vps s (kms -1 ) 


19.00±0.18 


31.8±0.7 


Age (Myr) 


57±5 


4.5±2.5* 



★Data from lTorres et al.H2010) 
fData fromglbrecht et a l1|2009| 
tData from Claret eTlu^feoiO} 

Note. — v synj denotes the stellar rotation speed for an aligned star which 
rotation is synchronized with the orbital period. vpSj denotes the pseudo 
synchronized value as de fined by [Hut (1981 1. To calculate vpSj we used a 
ftps of 0.8 (see |Hut|198f] Figure 37""^ 



5. THE ALIGNMENT IN CONTEXT 

Having established the absolute dimensions, age, and state 
of rotation in the EP Cru system we can now compare EP Cru 
to its apparently younger sibling DI Her. In Table [5] we reprint 
some of the values from EP Cru. According to these values 
the two systems are similar, apart from two characteristics: 
1) Their ages, EP Cru is about an order of magnitude older, 2) 
EP Cru appears to have aligned axes which is definitely not the 
case for DI Her. We would like to find a picture in which the 
misalignment in the young DI Her system can be explained as 
well as the alignment in the older EP Cru system. 

Knowing that the scaled radii are large enough in these 
systems to allow for substantial tides, we might suspect that 
the difference in spin-orbit alignment is a result of observing 
these systems at different stages in their evolution rather than 
them having two different formation and evolution paths. The 
hypothesis would be that both stars had misaligned axes, and 
we see EP Cru with aligned axes only because it is older and 
tides have had enough time to align the axes. 

The large eccentricities seen in both DI Her and EP Cru is 
consistent with this hypothesis, because tides first align and 
synchronize rotation and only on a longer timescale do they 
circularize the orbit. This is mainly due to the higher amount 
of angular momentum stored in the orbital motion compared 
to the stellar rotation, and for systems with a low-mass sec- 
ondary this is not necessary the case. However another finding 
makes the hypothesis difficult to reconcile with current tidal 
theories. The stars rotate at ~ 9 times the speeds expected for 
synchronized or pseudosynchronized states (Table [5]). Thus 
tides have not yet synchronized the stellar rotation speeds in 
the EPCru system. Formulations of tidal interactions pre- 
dict that damping of any significant spin-orbit misalignment 
should occur on the same time scale as synchronization of the 



EP Crucis 
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50 100 150 200 
rotational velocity (km s" 1 ) 

Fig. 7. — Tidal evolution of a system similar to EPCru, but with mis- 
aligned spin axes. The blue and red lines show the evolutions of primary 
and secondary obliquities (angle between stellar spin and orbital plane) in 
a system with the parameters of EP Cm. However we started the run with 
obliquities which have been measured in the DI Her system and faster stellar 
rotation. We included the stellar evolution (in particular the change of the 
stellar radii with time) of the system as estim ated with the Y 2 -evolutionary 
tracks ( Yi et al. 2001 Demarque et al. 2004) and set the viscous time (ty) to 
50 000 years, about 1 0UU times larger than what is normally assumed for late 
type stars. (A lower value of ty would lead to an overall faster tidal evolution 
but will leave the ratio of the alignment and synchronization timescales un- 
changed.) There is little evolution in the stellar obliquities until the rotation 
speeds approach the pseudo-synchronized value for rotation (Vps), which is 
indicated by the vertical line and is currently similar for the two stars (Ta- 
ble[5). 



rotation ( |Hut||l98l[ |Eggleton & Kiseleva-Eggleton|[200l| pl 
This is because in these tidal models, a single coefficient de^ 
scribes the coupling between tides and rotation. When stel- 
lar rotation is much faster than the synchronized value rota- 
tion around any axis is damped by about the same amount. 
Thus the angle between the overall angular momentum and 
stellar spin does not change: only the rotation speed is re- 
duced. When the stellar rotation around a axis parallel to the 
orbital angular momentum approaches the synchronized value 
than rotation around this axis couples less to tides. Rotation 
around any other axis is still damped by tides, which only 
ceases when the rotation around these axes stops. The stellar 
rotation aligns to the orbital axis. 

To illustrate this point we used the TOPPLE tidal-evo lution 
code developed by Eggleton & Kiseleva-Eggleton (2001). For 
this simulation we used the EP Cm parameters from Table B] 
but with initial obliquities taken from DIHer, and an initial 
faster stellar rotation speed at zero-age main sequence. The 
results are shown in Figure [7] The stellar obliquities remain 
large until the stellar rotation speeds approach synchroniza- 
tion, at which point obliquities are damped. This suggests that 



EP Cm had aligned axes when it was as young as DI Her, im- 
plying in turn that DI Her and EP Cm do not represent differ- 
ent stages of one evolution, but rather two different evolution 
paths. 

At the moment it is not possible to make more general state- 
ments as only a few measurements of obliquities have been 
carried out in close double-star systems. Furthermore most of 
these have been conducted in Algol systems which have un- 
dergone mass transfer (see Table 1 of JAlbrecht et al.||2011| ). 
Obliquity observations should be carried out in a variety of 
systems. Of particular interest would be young systems with 
short orbital periods with and without a third star. The sys- 
tems should be young to minimize the influence of tides, they 
should have orbital periods ranging from few days to few tens 
of days. Obliquity measurements in these systems would be 
helpful in testing predictions of KCTF and thereby of close 
binary formation. Measurements of obliquities in wider sys- 
tems would probe the length scale over which the primordial 
angular momentum was influential. Conducting such mea- 
surements is the aim of the BANANA project. 

6. SUMMARY 

We have analyzed high resolution spectra of the eclips- 
ing close double star system EP Cm. We obtained absolute 
dimensions and showed that the rotation axes of both stars 
are aligned with each other and the orbital rotation (/3 p = 
-1.8 ±1.6° and |/3 S | < 17°). 

EP Cm is similar in its orbital and stellar characteristics to 
DI Her. The two exceptions are that DI Her is younger and has 
two strongly misaligned stellar rotation axes. We have been 
unable to show that both systems represent different stages of 
one evolution path. This is because the stars in EPCru ro- 
tate at a few times their synchronized value and tidal theory 
predicts that synchronization occurs around the same time as 
alignment. Therefore the two systems likely represent two 
different formation scenarios rather then two different evo- 
lutionary stages. The sample of close double star systems 
for which the obliquities are measured remains small. We 
plan to ratify this situation by measuring obliquities in more 
close double star systems in the framework of the BANANA 
project. 
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